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Abstract 

The analytical solution of the Schrodinger equation for the Manning-Rosen potential plus a ring- 
shaped like potential is obtained by applying the Nikiforov-Uvarov method by using the improved 
approximation scheme to the centrifugal potential for arbitrary / states. The energy levels are 
worked out and the corresponding normalized eigenfunctions are obtained in terms of orthogonal 
polynomials for arbitrary / states. 
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I. INTRODUCTION 



As known, one of the main objectives in theoretical physics since the early years of 
quantum mechanics (QM) is to obtain an exact solution of the Schrodinger equation (SE) 
for some special potentials of physical interest. Since the wave function contains all necessary 
information for full description of a quantum system, an analytical solution of the SE is of 
high importance in non-relativistic and relativistic quantum mechanics l|, |2|. There are few 
potentials for which the SE can be solved explicitly for all n and / quantum states. 

The non-central potentials are needed to obtain better results than central potentials 
about the dynamical properties of the molecular structures and interactions. Researchers 
added ring shapeclpotentials to certain potentials, i.e. Coulomb, Kratzer [jsl] and Manning- 
Rosen potentials 41] to obtain non-central potentials. Ring-shaped potentials can be used 
in quantum chemistry to describe the ring shaped organic molecules such as benzene and 
in nuclear physics to investigate the interaction between deformed pair of nucleus and spin 
orbit coupling for the motion of the particle in the potential fields. 

It would be interesting and important to solve the SE for the Manning-Rosen potential 
plus a Ring-Shaped like potential for / 7^ 0, since it has been extensively used to describe 
the bound and continuum states of the interacting systems. Thus, one can obtain the 
energy eigenvalues and corresponding eigenfunctions of the one particle problem within this 
potential. The central Manning-Rosen potential is defined by 



a{a — l)exp{—2r/b) Aexp{—r/b) 



k = 2fi/ff. 



[1 — exp{—r /b)y {l — exp{—r/b) 
where A and a are dimensionless parameters, and b is the screening parameter. This po- 
tential is used as a mathematical modeling of the diatomic molecular vibrations and it 
constitutes a convenient model for other physical situations. It is known that, using for this 
potential the SE can be solved exactly for s-wave (/ = 0) jsj. Unfortunately, for an arbitrary 
/-states (/ ^ 0), the SE does not admit an exact solution. In such a case, the SE can be 

n 

solved numerically or approximately using approximation schemes [6[. 

Thepotential which is solved in this study is obtained by adding ring-shaped potential 
term j?] as. 



a{a — l)exp{—2r/b) Aexp{—r/b) 



+ 



f3cos9 

6^(1 — exp{—r/b)Y — exp{—r/b) ' ' 



+ 



:i.2) 
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So far, many methods were developed and app 



factorization 



ied, such as supersymmetry (SUSY) 



13, 



10| , Laplace transform approach llj and the path integral method [l2|, to 



solve the radial SE exactly or quasi-exactly for / ^ within these potentials. An other 
method known as the Nikiforov-Uvarov (NU) method [j] was proposed for solving the SE 
analytically. Many works show the power and simplicity of NU method in solving central 
and noncentral potentials 14j4l6|. This method is based on solving the second order linear 



differential equation by reducing to a generalized equation of hypergeometric type which is 
a second order homogeneous differential equation with polynomials coefficients of degree not 
exceeding the corresponding order of differentiation. 

In this study, we obtain the energy eigenvalues and corresponding eigenfunctions for 
arbitrary I states by solving the SE for the Manning-Rosen potential plus a ring-shaped like 
potential using NU method. Moreover, by chancing parameters we also obtain solutions 
for Hulten potential [l7|, and Hulten plus ring shaped potential. The NU method was also 
used to solve the SE for this pot ential. But our obtained results are more general and more 
efficient those given in Ref. [l8| 

The organization of this paper is as follows. The SE within Manning Rosen plus a ring- 
shaped like potential is provided in Section |Tll Bound state solution of the radial SE by 
NU method is presented in Section IIIII The solution of angle- dependent part of the SE is 
presented in Section IIVI and, the numerical results for energy levels and the corresponding 
normalized eigenfunctions are presented in Section |Vl Finally, some concluding remarks are 
stated in Section |VT1 



II. THE SCHRODINGER EQUATION WITH THE MANNING- ROSEN POTEN- 
TIAL PLUS A RING-SHAPED LIKE POTENTIAL 

The Schrodinger equation in spherical coordinates is given as 



Considering this equation, the total wave function is written as 



(2.1) 



V'(r,^,0) = i?(r)e(^^)$(0), 



(2.2) 
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where the polar angle solution is given by 

$(0) = i-e 



^"^■^ m = 0,±1,±2, 



(2.3) 



Thus for radial and azimuthal SE for Manning Rosen plus ring-shaped like potential are 

i?(r) = 0, (2.4) 



R"(r) + -R'(r) 
r 



2u ^ 1 Ae-'"/^ 1 a(a - l)e-^''/'' A 



0"(^) + cotee'ie) + 



l3' + l3cose\ 



m 



sin?6 



0(9) = 



(2.5) 



respectively. 



III. BOUND STATE SOLUTION OF THE RADIAL SCHRODINGER EQUA- 
TION. 



As know, Eq.(2.4) is the radial SE for Manning- Rosen plus a ring-shaped potential. In 
order to solve Eq.(2.4) with A = /(/ + 1) 7^ 0, we must make an approximation for the 
centrifugal term. When r/b << 1, we use an improved approximation scheme 19| to deal 
with the centrifugal term 



-r/b 



(1 _ e-r/b)2 



1 1 
12V 



1 

62 



-r/b 



(1 _ e-r/b)2 

(3.1) 



where the parameter Co = ^ is a dimensionless constant. However, when Co = the 
approximation scheme becomes the convectional approximation scheme suggested by Greene 
and Aldrich j20|. 

We assume R{r) = ^x(r) in Eq.(2.4) and the radial SE becomes 



X"{r) + 



2/i 1 Ae-^/^ _ 1 a(a - l)e-2^/^ _ A 
e'-l^) (1 - e-'-/'')2 " 62 



Co + TT—Tn^. X(r) = 0. 

(3.2) 

In order to transform Eq.(3.2) in to the equation of the generalized hypergeometric-type 



which is in the form 



13| 



r 



x\s) + -x'(s) + ^x(s) = 



we use the transformation s = e~'~/^. Hence we obtain 



(3.3) 
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s{l-s) 



-e^(l - sf + As{l - s) - a{a - l)s^ - 

s 



{1-sy 



X{s) = 0, (3.4) 



where we use the foUowing notation for bound states 



— e 



2// 



Eh'^, E <0. 



(3.5) 



Now, we can successfully apply NU method of definition for eigenvalues of energy. By 
comparing Eq.(3.4) with Eq.(3.3) we can define the foUowings 

f(s) = 1 — s, a{s) = s(l — s), 

a{s) = s^f-e^ -A-a{a-l)- Acq] + s[2e^ + A + 2\Cq - A] + [-e^ - Acq]. (3.6) 



If we take the following factorization 



X{s) = 0(s)y(s), 



(3.7) 



for the appropriate function <p[s) the Eq.(3.4) takes the form of the well known 
hypergeometric-type equation. The appropriate (j){s) function must satisfy the following 
condition 

As) Ms) 



0(s) a{s) 



where function 7r(s) is defined as 



(3.9) 



Finally the equation, where y{s) is one of its solutions, takes the form known as 
hypergeometric-type 

a{s)y"{s) + T{s)y'{s) + Xy{s) = 0, (3.10) 



where 



and 



X = k + TT , 



t{s) = f(s) + 27r(s). 



(3.11) 



(3.12) 
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For our problem, the 7r(s) function is written as 



^[s) = — ± ^/s'^[a - k] - s[b - k] + c, 



(3.13) 



where the values of the parameters are 

1 



a=| + e^ + A + a{a - 1) + ACq, 
6 = + A + 2ACo - A, 



c = + ACn 



The constant parameter k can be found complying with the condition that the discrimi- 
nant of the expression under the square root is equal to zero. Hence we obtain 



A;i 2 = (6 - 2c) ± 2Vc2 + c(a-6). 
Now we can find four possible functions for 7r(s) 



(3.14) 



, , -s \ (V^- Vc + a-6)s-x/5 /or A; - (6 - 2c) + 2 ^c^ + c(a - 6), 

7r(s) = — ± < (3.15) 

2 I (\/c + \/c + a-6)s -^c for k = {b - 2c) - 2v/c2 + c(a - 6). 

According to NU method, from the four possible forms of the polynomial 7r(s), we select 
the one for which the function t(s) has the negative derivative. Therefore, the appropriate 
function 7r(s) and t(s) are 



7r(s 



) = Vc- s 



- + \/c+ Vc + a-b 



t{s) = 1 + 2^/c - 2s 1 + Vc + a - &] , 



for 



A; = (6 - 2c) - 2^yc'^ + c{a-b). 
Also by Eq.(3.11) we can define the constant A as 



(3.16) 
(3.17) 
(3.18) 



A = 6 - 2c - 2y/c'^ + c{a-b) 



- + -\/c+ Vc + a - b 

Zi 



(3.19) 



Given a nonnegative integer n, the hypergeometric-type equation has a unique polyno- 
mials solution of degree n if and only if 
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A = A„ = —riT 



a ,(ra = 0,l,2...), 



(3.20) 



and Am 7^ A„ for m = 0, 1, 2, n — 1 211], then it follows that 



\n, = h-2c- + c(a - h) - 



- + a/c+ Vc + a - 6 



2?T,r- 



1 + ^a/c + a/c + a — 6 j + nr(rar — 1 



(3.21) 



We can solve Eq.(3.21) explicitly for c and by using the relation c = e^ + ACo which brings 



-i 2 



A + 1/2 + A(l + 2?2r) + n,.(?2r + 1) - A 
2A + 1 + 2n 

where A = a/1/4 + a(a - 1) + A. 

We substitute into Eq.(3.5) with A = /(/ + 1) which identifies 



ACn 



(3.22) 



Enr,l 



Tlr + 1/2 + 



(/ - nr){l + Tlr + I) - A 



l{l + l)Co 



(3.23) 



2A + 1 + 2nr 

Now, using NU method we can obtain the radial eigenf unctions. After substituting 7r(s) 
and a{s) into Eq.(3.8) and solving first order differential equation, it is easy to obtain 



s = s 



(3.24) 



where = 1/2 + A. 

Furthermore, the other part of the wave function y(s) is the hypergeometric-type function 
whose polynomial solutions are given by Rodrigues relation 



yn[S) 



(3.25) 



p(s) ds^ 

where Bn is a normalizing constant and p{s) is the weight function which is the solutions of 
the Pearson differential equation. The Pearson differential equation and p{s) for our problem 
is given as 



(ap)' = TP, 



(3.26) 



Pis) 



(3.27) 
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respectively. 

Substituting Eq.(3.27) in Eq.(3.25) we get 



; 2k—l+nr 



Then by using the following definition of the Jacobi polynomials 22 1 



-1)" 



n!2"(l - + sfds' 



[{1 - + 



we can write 



and 



[s-+-^l - 5)''+"] = Cns^il - sfPt''\l - 2s). 
ds^ '- 

If we use the last equality in Eq.(3.28), we can write 

l/„.(s)=C„p(^v/^'2^-i)(l-2.). 
Substituting p(s) and UnXs) into Eq.(3.7), we obtain 



X 



Using the following definition of the Jacobi polynomials 22 1 



P^''f'\s) = F i -n,a + b + n + 1,1 + a; 



n!r(a + l) 21 V ■ 2 

we are able to write Eq.(3.33) in terms of hypergeometric polynomials as 



rirli \2^c + 1) 21 ^ ' 
The normalization constant can be found from normalization condition 



OO 1 

|i?(r)|Vrfr = y \x{r)\^dr = b J ^\x{s)\^ds = 1, 
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by using the following integral formula 



23| 



;i - ^)2(<5+l)22A-l ) ^(_^^^ 2(5 + A + 1) + rir, 2A + l]z)\ dz 



21 



_ {nr + 5 + l)n,!rK + 2(5 + 2)r(2A)r(2A + 1) 
~ K + 5 + A + l)rK + 2A + l)r(2(5 + A + 1) + n,) ' ^ ' ^ 

for (5 > ^ and A > 0. After simple calculations, we obtain normalization constant as 



^ ^ nr\2^c{nr + K + v/^)r(2(ir + y/^) + n,) 

^' 6K + ir)rK + 2v^+ i)r(n, + 2K) • ^ ' 

IV. SOLUTION OF AZIMUTHAL ANGLE-DEPENDENT PART OF THE 
SCHRODINGER EQUATION 

We may also derive the eigenvalues and eigenvectors of the azimuthal angle dependent 
part of the SE in Eq.(2.5) by using NU method. The boundary condition for Eq.(2.5), O(^) 
require to be taken as a finite value. Introducing a new variable x = cos6, Eq.(2.5) is brought 
to the form 

e"(x) - -^e'(x) + — [A(l - x') - - (/?' + e(x) = 0. (4.1) 

1 — a;^ (1 — a;^)^ 

After the comparison of Eq.(4.1) with Eq.(3.3) we have 

f{x) = -2x ,a{x) = l- x'^ ,a{x) = -Xx^ - f3x + {X-m^ - /3'). (4.2) 
In the NU method the new function 7i{x) is calculated for angle- dependent part as 

tt{x) = ±^/x^{X -k) + (5x-{X- [5' -m^ - k). (4.3) 

The constant parameter k can be determined as 

2X-m^-P' u 
= ± -, (4.4) 

where u = a/ {m'^ + /3')^ ~ P"^- 

The appropriate function 7r(x) and parameter k are 



tt{x) 



m? + /?' + M jm? + l3' — u 
x\l t: h 



(4.5) 
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k 



2A - - u 



(4.6) 



2 2 

The following track in this selection is to achieve the condition r' < . Therefore t{x) 
becomes 



t{x) — —2x 



m? + /3' + -u 



m? + (5' — u 



(4.7) 



We can also write the values A = A; + 7r'(s) as 



2A - - w? 



u 



rn? + P' + u 



(4.8) 



also using Eq.(3.20) we can equate 



A 



2A - ^' - u m? + /3' + m 



N 



2N 



1 + 



m? + /3' + u 



2 2 V 2 

In order to obtain unknown A we can solve Eq.(4.9) exphcitly for A = /(/ + 1) 



+ N{N-1). (4.9) 



(4.10) 



where C = y^^^^, and 



A = + C + 2iVC + iV(iV + 1) = (AT + C)(iV + C + 1) = /(/ + 1), (4.11) 



then 



l^N + C- (4.12) 

Substitution of this result in Eq.(3.23) yields the desired energy spectrum, in terms of rir 
and I quantum numbers. Similarly, the wave function of polar angle dependent part of SE 
can be formally derived by a process to the derivation of radial part of SE. 



0(x) = (l-a;)(^+^)/2^ 



(4.13) 



p(x) = (l-x)^+^(l + a;) 



\B-C 



(4.14) 



d 



y^{x) = B^{1 - x)-(^+^)(l + xf-'^^ [{1 - x)^+^+^(l + x)^-^+^] , (4.15) 
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where 



^ jm? + [5' + u ^ I m? + /3' — u 



2 ' V 2 
From the definition of Jacobi polynomials, we can write 



^ + = (-l)^2^(l-x)^+^(l + x)^-^Pi^+^'^-^)(a;). (4.16) 



Substitution of Eq.(4.16) into Eq.(4.15) and after long but straightforward calculations we 
obtain the following result 

e^(x) = C^(l - x)(^+^)/2^1 + x)(^-^)/V;^+^'^-^)(x), (4.17) 

where Cm is the normalization constant. Using orthogonality relation of the Jacobi polyno- 



mials 



22 1 the normalization constant can be found as 



^ , i2N + 2B + l)T{N + mN + 2B + l) 

\l 2^B+iY(^N + B + C+l)T{N + B-C + l)' ^ ' 

V. NUMERICAL RESULTS AND DISCUSSION 

Solution of the SE for the Manning-Rosen potential plus a ring-shaped like potential are 
obtained by applying the Nikiforov-Uvarov method in which we used the improved approxi- 
mation scheme to the centrifugal potential for arbitrary / states. The energy eigenvalues and 
corresponding eigenfunctions are obtained for arbitrary / quantum numbers. Two important 
cases must be emphasized in the results of this study. In the first case which /3 = /?' = 
the potentials turn to central MR potential. For this case, by using u = m? , ( = \m\ and 
/ = + \m\ {N = 0, 1, 2...) then I > \m\ by substituting this / values in Eq.(3.23) we obtain 
energy spectrum for MR potential. These results are consistent with those of works in Zhao- 



You Chen and et al, 



24( 1 . Also, if a = or a = 1, then in this case, for 6 = 1/b, A = 2b 



potential turns to the Hulten potential with / > |m|. This energy values are consistent 
with those of works in Varshni The second case is the general situation where P ^ 
or /3' 7^ 0. By changing /3 and /3' we obtain energy values of MR plus different type ring 
shaped like potentials. In Table l(for a = 0.75) and Table 2(for a = 1) we show energies of 
the bound states for the Manning-Rosen potential plus ring shaped potentials for different 
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approximation to it ^ = [Cq + ^^lj-rihyi \ (dashed) and the conventional approximation to 
it ^ = 73- [tt^t^t/fw] (long dashed) as function of the variable r are displayed from figures 1 to 



values of ^, m, Z, n and 1/h — 0.025, A — 2b, where n — rir + I + 1, is the usual principle 
quantum number and rir is the number of nodes of the radial wave functions. To show 
accuracy of our results, the plots of the centrifugal term (solid line), the improved new 

^ [Co + (YTF^Vfcja] 

b2'[(l_e-r/b)2j 

3 with different potential range parameter b=l,2,4. It is shown that the our approximation 
Eq.(3.1) is a good approximation to the centrifugal term for short range potential range, i.e. 
large b. 

Finally we want to deal with some restrictions about bound state solutions of SE for 

Manning Rosen plus ring shaped like potential. Firstly, it is seen from Eq(4.4) and expression 
from u that in order to obtain real energy values the condition (m^ + 2/3')^ > 4/3^ must be 
hold. Since the parameters (3 and f3' are real and positive, we can write 

> 2(/3-/3'). (5.1) 

li /3 < /3' the inequality in Eq.(5.1) is provided automatically. But if ^ > ^' then m 
becomes bounded. Secondly, in Eq.(3.23) if 



l{l + l)Co < 



71^ + 1/2 + 



(/ - nr){l + Hr + 1) - A' 



2 



(5.2) 



2A + 1 + 2nr 

then energy eigenvalues take nonnegative values, this means there is no bound states. If we 

take Co — for the approaches in centrifugal term as some previous studies, the restriction 

on quantum numbers is removed. Finally from Eq (3.21) we obtain 

^ X+1/2 + A{l + 2nr) + nr{nr + l) - A 
^ = -(2A + l + 2n.) ■ ^^-^^ 

For bound states since c > (so > ) and using the fact that 2A + 1 + 2nr > 0, we 

obtain 

^ - 1/2 - /(/ + 1) - A > riririr + 2A + 1), 

A>^ + l{l + l)+A. (5.4) 

If both conditions in Eqs.(5.1-5.2 and 5.4) are satisfied simultaneously, the bound states 
exist. Thus the energy spectrum equation in Eq.(3.23) as limited, i.e, we have only the 
finite numbers of energy eigenvalues. Figure 4 and Figure 5 shows the region of possible 
values of rir a-nd I quantum numbers according to our analysis in Eq.(5.2) and Eq.(5.4) for 
A = 80, q; = 1 and a — 0.75 respectively. 
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VI. CONCLUSION 



Analytical calculations of energy eigenvalues for an arbitrary / state and corresponding 
eigenf unctions in the Manning-Rosen potential plus a Ring-Shaped like potential is done 
by using Nikiforov-Uvarov Method in this paper. The energy eigenvalue expression for 
Manning- Rosen potential plus a Ring-Shaped like potential is given by Eq.(3.23). By some 
specific value of a, /3 and /3' parameters one can observe that the some results of previous 
studies for finding MR and Hulten potential can be found. We also obtain some important 
restrictions on quantum numbers about bound state solutions of SE. We can conclude that 
our results are not only interesting for pure theoretical physicist but also for experimental 
physicist, because the results are exact and more general. We have also examined that, 
obtained results are different from results given in 18 1. Therefore, calculation in ji~ 



should be checked and recalculated and compared with our results and 17|, |24 |. 
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3.119190 


-0.030831 


1 


1 


1 


1 





1.618034 


3.618034 


-0.027392 


1 


1 


1 





1 


2.618034 


3.618034 


-0.026949 





1 


1 


1 





1.732051 


3.732051 


-0.025026 



/3 


/3' 


m 


N 




I 


n 


E 





1 


1 





1 


2.732051 


3.732051 


-0.024630 











3 








4 


-0.024017 


1 





2 


1 





1.931852 


3.931852 


-0.021403 


1 





2 





1 


2.931852 


3.931852 


-0.021065 








1 


2 





1 


4 


-0.020809 


1 


1 





2 





1 


4 


-0.020809 











1 


2 


2 


4 


-0.020299 


1 


1 





1 


1 


2 


4 


-0.020299 


1 





3 








2.981188 


3.981188 


-0.020271 








1 





2 


3 


4 


-0.019976 


1 


1 








2 


3 


4 


-0.019976 


1 


1 


3 








3.302776 


4.302776 


-0.015781 





1 


3 








3.316625 


4.316625 


-0.015611 





1 





2 





1.414214 


4.414214 


-0.015053 





1 





1 


1 


2.414214 


4.414214 


-0.014733 


1 


1 


2 


1 





2.414214 


4.414214 


-0.014733 


1 


1 


2 





1 


3.414214 


4.414214 


-0.014463 





1 








2 


3.414214 


4.414214 


-0.014463 





1 


2 


1 





2.449490 


4.449490 


-0.014336 



TABLE I: Energies of the bound states for the Manning Rosen potentials plus ring shaped potentials 
for different values of /?, /?', m, i, and a = 0.75,1/6 = 0.025, A = 26 calculated using Eq.(4.12) 
and Eq.(3.23) for = /i = 1. 
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p 


/3' 


m 


N 




I 


n 


E 




















1 


-0.487578 











1 








2 


-0.112812 








1 








1 


2 


-0.112760 


1 


1 











1 


2 


-0.112760 





1 











1.414214 


2.414214 


-0.073653 


1 


1 


1 








1.618034 


2.618034 


-0.060874 





1 


1 








1.732051 


2.732051 


-0.054947 


1 





2 








1.931852 


2.931852 


-0.046192 











2 








3 


-0.043759 








1 


1 





1 


3 


-0.043707 


1 


1 





1 





1 


3 


-0.043707 








2 





2 


2 


3 


-0.043602 


1 


1 








1 


2 


3 


-0.043602 





1 





1 





1.414214 


3.414214 


-0.031215 


1 


1 


2 








2.414214 


3.414214 


-0.031089 





1 








1 


2.414214 


3.414214 


-0.031089 





1 


2 








2.419490 


3.119190 


-0.030230 


1 


1 


1 


1 





1.618034 


3.618034 


-0.026609 


1 


1 


1 





1 


2.618034 


3.618034 


-0.026473 





1 


1 


1 





1.732051 


3.732051 


-0.024363 



/3 


/3' 


m 


N 




/ 


n 


E 





1 


1 





1 


2.732051 


3.732051 


-0.024221 


1 





2 


1 





1.931852 


3.931852 


-0.020903 


1 





2 





1 


2.931852 


3.931852 


-0.020750 











3 








4 


-0.020000 


1 





3 








2.981188 


3.981188 


-0.019975 








1 


2 





1 


4 


-0.019948 


1 


1 





2 





1 


4 


-0.019948 


1 


1 





1 


1 


2 


4 


-0.019844 








1 


1 


1 


2 


4 


-0.019844 








3 








3 


4 


-0.019687 


1 


1 








2 


3 


4 


-0.019687 


1 


1 


3 








3.302776 


4.302776 


-0.015583 





1 


3 








3.316625 


4.316625 


-0.015417 





1 





2 





1.414214 


4.414214 


-0.014594 





1 





1 


1 


2.414214 


4.414214 


-0.014468 


1 


1 


2 


1 





2.414214 


4.414214 


-0.014468 


1 


1 


2 





1 


3.414214 


4.414214 


-0.014290 





1 








2 


3.414214 


4.414214 


-0.014290 





1 


2 


1 





2.449490 


4.449490 


-0.014082 



TABLE II: Energies of the bound states for the Manning Rosen potentials plus ring shaped poten- 
tials for different values of (3, 13' ,m, I, rir and a = 1,1/b = 0.025, A = 2b calculated using Eq.(4.12) 
and Eq.(3.23) foi ji = h = 1. 
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2 4 6 8 10 

r 

FIG. 1: The plots of the centrifugal term (solid line) the improved new approximation to it (dashed) 
and the conventional approximation to it (long dashed) as the function of the variable r with 
potential range parameter b=l 




b=2 
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FIG. 2: The plots of the centrifugal term (solid line) the improved new approximation to it (dashed) 
and the conventional approximation to it (long dashed) as the function of the variable r with 
potential range parameter b=2 
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0,0 




FIG. 3: The plots of the centrifugal term (solid line) the improved new approximation to it (dashed) 
and the conventional approximation to it (long dashed) as the function of the variable r with 
potential range parameter b=4 
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FIG. 5: The region of possible values of and I for a = 1, l/b=0.025 and A=2b 
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